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SECTION  1 
INTRODUCTION 

Ihert'  is  evidisnce  that  the  boundary  layer  on  a reentry  nosetip  can  be  turbu- 
lent in  the  neighborhood  of  the  sonic  point  and  laminar  or  ''laininar-1  ike"  in  a down- 
stream reyion  near  the  sphere-cone  tangent  point.  This  phenomenon  is  commonly  called 
relaniinari^ation  or  laminarization  and  re^.'lts  in  locally  reduced  heat  transfer  rates. 
The  major  impact  of  this  locally  reduced  heating  is  expected  to  be  on  the  thermostruc- 
tural  performance  of  a nosetip.  In  order  to  evaluate  the  impact  of  this  effect  on 
nosetip  performance,  it  is  necessary  to  develop  a method  for  predicting  laniinarization 
and  evaluating  the  associated  heat  transfer  distributions  that  is  compatible  with  the 
codes  used  for  r.usetip  design  and  performance  studies.  This  methodology  and  a pre- 
liminary sensitivity  study  are  the  primary  objectives  of  the  present  study. 

Reiaminarization  has  been  observed  in  strongly  accelerated  axisyinmetric  and 
two-dimensional  flows  as  well  as  in  flows  with  convex  streamwise  curvature.  Signifi- 
cant relaminarization  and  reductions  in  local  heat  trans'er  rate  have  been  inferred 
from  measured  temperature  distributions  on  calorimeter  models  tested  in  the  Naval 
Ordnance  Laboratory  Tunnel  No.  fl  and  on  ATJ-S  graphite  nosetip  models  tested  in  the 
Air  Force  Materials  Laboratory  50  MW  arc.  Since  the  phenomenon  exists  for  such  a 
wide  range  of  conditions  and  configurations,  it  is  necessary  to  delineate  the  parame- 
ters most  relevant  to  nosetip  laminarization.  This  cnoicc  is  based  on  an  analysis 
of  the  mechanics  of  nosetif)  boundary  layers  and  on  the  need  to  develop  a prediction 
scheme  suitable  fur  inclusion  in  the  shape  change  codes  used  for  design  and  perfor- 
mance studies. 

Whereas  it  is  anticipated  that  the  major  impact  of  relaminarization  on  nose- 
tip  perfopiiance  is  thermostructural , the  scope  of  the  present  investigation  is  limited 
tc  the  evaluation  of  the  associated  heat  transfer  distribution.  This  objective  has 
been  accomplished  by  the  following  steps 

a.  Review  and  evaluation  of  theoretical  and  experimental  work  on  laminariza- 


b.  Development  and  incorporation  into  a design  code  of  a relaminarization 
heat  transfer  correlation 


Preceding  page  blank 
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c.  Sensitivity  study  to  determine  the  effects  of  trojectory  and  nose  radius 
variation. 

The  results  of  the  sensitivity  study  indicate  significantly  reduced  heating 
rates  near  the  sphere-cone  tangent  point  of  a nosetip  and  indicate  a need  for  further 
work  in  order  to  determine  the  full  extent  of  the  thermostructural  ramifications. 

Section  2 of  this  report  highlights  some  of  the  important  aspects  of  the  lit- 
erature review  and  evaluation  published  separately  (Wolf  and  Anderson,  Reference  1), 
Section  3 documents  the  methods  and  results  of  the  correlation  development  task,  and 
Section  4 details  the  results  and  conclusions  of  the  sensitivity  study.  Primary 
conclusions  are  presented  in  Section  1. 


SECIION  2 
LITERATURE  SURVEY 

Ttie  objective  of  this  ph.ise  of  the  jnvestiqation  was  to  assemble  and  evaluate 
available  information  on  the  phenomenon  of  relaminar  ization.  Ttie  detailed  results 
of  this  critical  survey  are  reported  by  Wolf  and  Anderson  (Reference  1). 

Since  what  were  probably  the  first  observations  of  laminarization  by  Sternberg 
in  1954  (Reference  2),  a large  number  of  theoretical  and  experimental  studies  have 
been  devoted  to  the  subject.  These  references  were  evaluated  in  terms  of  the  need  to 
implement  the  results  in  the  shape  change  codes  used  for  nosetip  desigti.  Hence,  the 
choice  of  variables  is  limited  to  the  boundary  layer  parameters  that  appear  in  the 
von  Kan.ian  integral  momentum  equation. 

Known  lamina-izing  flows  can  be  broadly  classified  into  two  groups,  which 

are; 

a.  Boundary  layers  subjected  to  large  streamwise  accelerations 

b.  Boundary  layers  subjected  to  large  normal  accelerations. 

Of  these  two  groups,  (a)  has  been  the  most  widely  studied.  Experimental  data  covering 
a wide  range  of  pressure  gradients,  initial  Reynolds  numbers,  configurations  and  heat 
and  mass  transfer  rates  are  available.  Some  of  the  data  are  very  detailed,  including 
measurements  of  the  structure  of  the  turbulence  and  the  inlormi ttency  of  a laminari- 
zing  boundary  layer.  In  contrast,  there  is  relatively  less  data  on  boundary  layers 
v/hich  laminarize  as  a result  of  large  normal  accelerations.  In  this  case  there  are 
only  a few  reported  studies  and  the  experiments  are  less  detailed.  This  information 
is  summarized  on  Figure  4 of  Reference  1,  which  is  a tabulation  of  experimental  sour- 
ces and  the  parameters  that  were  measured. 

A similar  situation  exists  with  respect  to  theoretical  descriptions  of  lamin- 
arizing  flows.  The  theories,  (see  Wolf  and  Anderson,  Reference  1)  range  from  methods 
based  on  modifications  of  the  usual  approach  to  the  integral  momentum  equation  to 
multi-equation  turbulence  models.  Some  of  the  theories  include  compressibility  and 
streamwise  curvature  as  well  as  the  necessary  extensions  for  the  prediction  of  heat 
transfer.  Following  the  suLmittal  of  Reference  1,  several  recent  additional  contribu- 
tions by  harasimha  and  Sreenivasan  (Reference  3),  Baker  and  Launder  (References  4 and 
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5),  arid  Uoissler  (References  6 and  7)  were  examined.  However,  the  main  conclusion 
of  Wolf  and  Anderson  (Reference  I)  must  remain  as  stated,  namely,  that  there  is  no 
I'eally  reliable  and  well  tested  method  for  predicting  laminarizing  flows.  Even  if 
it  IS  jidged  that  one  (or  more)  of  the  existing  theories  shows  satisfactory  compari- 
son with  a certain  subset  uf  the  data,  it  is  extremely  unlikely  to  be  applicable  to 
a hypersonic  nosetip.  Furthermore,  the  kind  of  correlation  that  would  be  of  use  for 
practical  design  calc.ildtions  was  unavailable. 

A large  number  of  different  parameters  for  characterizing  relaminarization 
have  been  suggested.  Wolf  and  Anderson  (Reference  1)  evaluated  these  with  respect 
to  the  needs  of  the  present  study.  Their  conclusions  are  based  mostly  on  the  fact 
that  only  relutively  simple  integral  and  edge  variables  are  available  from  the  cur- 
rent design  codes,  but  also  on  the  lack  of  definitive  evidence  for  the  validity  of 
any  one  of  the  proposed  paraineter(s) . 


I 


I: 


i. 


SECTION  3 

CORRELATION  DEVELOPMENT 

It  is  inanediately  apparent  that  a nosetip  environment  con*^ains  both  of  the 
factors  known  to  cause  iaminarizatlon  of  other  flows,  i.e.,  streamwise  and  normal 
dccel ei-ation.  In  addition,  a nosetip  boundary  layer  is  influenced  by  the  complica- 
tions of  large  surface  roughness,  compressibility  and  surface  mass  transfer.  There- 
fore, it  is  necessary  to  identify  the  dominant  nosetip  laniinarization  mechanism  as 
a preliminary  step  in  developing  a correlation.  The  approach  to  this  task  is  to 

examine  separately  each  class  of  known  laminarizing  t'lows  and  on  the  basis  of  these 

results  construct  a correlation  for  a nosetip  boundary  layer. 

fhere  are  our  relatively  recent  ■ >rinienta1  studies  of  plane  boundary  layers 
with  streaimvise  acceleration  that  contai-  .itably  detailed  information  on  integral 
boundary  layer  parameters:  Rlackwelder  a.  r Kovasznay  (References),  Patel  and  Hoad 

(Reference  9),  Badri  Narayanan  and  Ramjee  (Reference  10),  and  Jones  and  Launder  (Ref- 
erence II).  The  data  from  these  experiments  were  examined  in  detail,  and  the  methods 

aivl  results  are  given  in  Section  3.1  below.  The  effort  to  correlate  these  data  is 

on'y  partially  successful  due  to  what  are  believed  to  be  inconsistencies  in  the  data. 
However,  sufficient  qualitative  information  was  obtainr..  to  justify  the  use  of  a cor- 
relation based  on  other  data,  as  described  in  Sectioti  3.2. 


3. 1 


DATA  REDUCTION  AND  CONSISTENCY  ANALYSIS 


in  order  to  facilitate  the  necessary  analysis,  the  data  wer-’  transferred  to 
(iiinrtied  cards  and  a code  was  written  to  test  various  aspects  nf  the  data  and  for  use 
as  a tool  in  correlation  development.  Input  to  the  code  are  the  measured  distribu- 
tions of  all  the  variables  in  the  integral  inomentuin  equation 


dRe  , 


Re  dS  ^ '^f/' 


i 


1 


where  s - streamwise  roordinate 

u^_  - boundary  layer  edge  velocity 

'*  = dis[)lacement  thickness 
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momentum  thickness 


T = wall  shear  stress 
w 


= edge  mass  density 


Vg  = edge  kinematic  viscosity 


Re  = Reynolds  number 


Re^j  = UgQ/Ug  = momentum  thickness  Reynolds  number 


H = 5*/0  = shape  factor 


K = \'g/Ug  dUg/ds  = pressure  gradient  parameter 


■ skin  friction  coefficient 


One  basic  section  in  the  code  generates  a cubic  spline  fit  to  any  variable  with 
a controllable  amount  of  smoothing  using  a least  squares  cr  terion.  The  method  is  an 
adaption  of  that  due  to  Reinsch  (Reference  12)  and  is  necessary  in  order  to  obtain 
accurate  interpolation  between  data  points  and  to  calculate  well  behaved  derivatives 
of  the  data  where  needed.  An  example  of  the  results  of  using  this  routine  is  shown 


in  Figure  1 where  a fit  to  the  measured  edge  velocity,  {Ug(s)),  is  compared  to  the 


data.  Figure  2 shows  the  corresponding  distribution  of  K(s),  calculated  from  the 
smoothed  cubic  spline  and  compared  to  that  given  by  Blackwelder  and  Kovasznay  (Ref- 
erence 8).  It  is  seen  that  the  method  gives  good  results,  both  in  terms  of  the  fit 
to  th-'  data  and  for  calculating  derivatives.  The  differences  between  the  curves  and 
the  data  indicated  by  Figures  1 and  2 are  not  significant  in  terms  of  the  effect  on 
the  subsequent  analysis  or  in  terras  of  evaluating  the  method.  Instead,  the  "goodness 
of  fit"  is  dependent  on  control  inputs  to  the  program,  v/hich,  when  appropriately  cho- 
sen, will  generate  cubic  splines  that  pass  exactly  through  the  data.  An  example,  of 
very  little  smoothing  is  given  in  Figure  3,  which  shows  the  shape  factor  distribution 
for  the  data  of  Blackwelder  and  Kovasznay  (Reference  8)  and  its  analytic  representa- 
tion. 


Another  section  of  the  code  contains  the  logic  and  subroutines  necessnM'  for 
using  the  data  to  solve  the  integral  momentum  equation.  For  example,  with  Kt,  . 

C^.  and  Re,^(0)  taken  from  the  data,  one  option  calculates  Re-j(s)  by  integration  of 
the  differential  equation  using  a fourth  order  Runge-Kutta  sr^-.me.  Another  option 
used  fits  to  the  measured  distributions  of  Re,  H and  Re,.(s)  ‘■q  calculate  C^(s).  The 
code  also  facilitates  the  prediction  of  boundary  layer  development  when,  as  is  the 
case  here,  the  initial  and  freestteam  conditions  are  known  and  when  a correlation 
that  specifies  a relationship  between  Re,^,  and  H is  postulated. 

The  results  of  the  analysis  of  the  data  and  various  correlations  are  given 
in  the  next  sections,  wherein  first  a special  case  and  then  the  general  case  of  stream- 
wise  accelerated  boundary  layers  are  considered. 
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Comparison  of  the  shape  factor  distribution  measured 
by  Biackwelder  and  Kovasznay  (Reference  8)  with  the 
smoothed  cubic  spline  fit. 


3.1.1  Et|ui  librium  Sink  i~lows 


An  equilibrium  sink  flow  is  a boundary  layer  for  which  the  pressure  gradient 


parameter, 


This  condition  is  obtained  experimentally  from  the  flow  between  convergent  plane 
walls.  Its  significance  is  that  is  describes  the  asymptotic  behavior  of  an  acceler- 
ated boundary  layer  (after  relaxation  effects)  subject  to  constant  K.  Therefore,  an 
equilibrium  sink  flow  represents  a limiting  state  for  an  accelerated  turbulent  boun- 
dary layer. 


For  this  circumstance,  it  has  been  shown  theoretically  and  experimentally  by 
Lainuer  and  Jones  (Reference  13),  and  Jones  and  Launder  (Reference  11)  that  K is  the 
single  independent  parameter  that  characterizes  the  flow.  A correlation  based  on 
these  studies  was  developed  and  tested  against  theory  and  experiment  with  satisfac- 
tory results.  This  correlation  is  given  by  the  following  expressions: 


Re,,^  = 0.0117 


Hf  = H^/1  . K/K, 


Turbulent  ; K £ Kj 


Re,,^  = 0.3762  K'-' 


H,  = 2.069 


Laminar  ; K £ K|^ 


Re.,  = (1  - W)  Re^^,  ■■  WRCf, 


H - (1  - W)  ^ Wll|^ 


(1  - cos  'iZ)/2 


' Laminar-Turbulent  ; ICj. 


< K K, 


log  K - log  K.J- 
To'g”K|^~  -’log’l'.^ 


The  parameters  chosen  for  this  correlation  are: 


= 0.68  (10-'’) 


= 8.0  do-'-) 


Kp  = 8.0  (lO'd 


1.2  to  1.5 


H 


fj 


where  is  the  experimental  shape  factor  at  the  upstream  position  where  K = 0. 


This  correlation  is  a composite  of  existing  theoretical  and  experimental  re- 


sults. When  K > K|^  it  gives  results  identical  to  the  known  exact  (analytical)  solu- 


tion for  a laminar  sink  flow  (see  Schlichting,  Reference  14)  For  K < Kj,  it  matches 
the  results  of  the  calculations,  based  on  a mixing  length  mode'',  of  Launder  and  Jones 
(Reference  13).  When  Kj  ^ K < the  given  form  of  the  weighting  factor,  W,  inter- 
polates between  the  turbulent  and  laminar  extremes  with  an  analytical  form  that 
matches  the  slopes  of  the  extremes  at  the  limit  points. 


When  K is  specified,  H and  Re^  are  determined  by  the  above  correlation,  and 
the  skin  friction  coefficient  is  found  from  the  integral  momentum  equation 


= 2(H  + 1)  RegK 


which  is  the  appropriate  form  when  dROg/dx  = 0. 

Results  obtained  from  the  correlation  and  representative  data  are  shown  in 
Figures  4,  5 and  6.  It  is  seen  that  the  trend  and  the  magnitude  of  the  data  are 
correctly  given  by  the  correlation. 


Therefore,  a simple  computational  scheme  has  been  established  for  determining 
the  asymptotic  state  (and  integral  properties)  of  an  initially  turbulent  boundary 
layer  when  the  imposed  pressure  gradient  is  such  as  to  make  K = constant.  This  cor- 
relation does  not  describe  the  process  of  evolution  to  the  asymptotic  state  but  it 
is  useful  for  estimating  the  magnitude  of  the  pressure  gradient  necessary  to  even- 
tually obtain  a given  degree  of  laminarization.  Thus,  Figure  4 shows  that  the  often 
quoted  value  of  K = 3(10'^)  corresponds  to  an  asymptotic  state  that  is  about  "half 
laminar"  and  "half  turbulent".  Furthermore,  the  summary  of  experimental  and  theoret- 
ical results  represented  by  the  correlation  shows  that  laminarization  is  a process 
that  occurs  to  varying  degrees  over  a range  of  pressure  gradients  given  by  0.7  (10"'’) 
K ''  8 (10-''). 


These  results  provide  useful  guidelines  for  identifying  many  qualitative  fei- 
tures  of  other  experimental  circumstances,  and  in  particular  are  good  indicators  of 
nonequilibrium  effects. 


3.1.2  Nonequi 1 ibrium  P lane  Flows 


A turbulent  boundary  layer  responds  to  changes  in  edge  conditions  by  processes 
that  require  finite  time.  When  conditions  are  such  that  a boundary  layer  is  continu- 
ously adjusting  to  changes  in  external  pressure  gradient,  it  is  nonequi 1 ibrium  if  the 
rate  of  adjustment  is  slower  than  the  rate  at  which  conditions  change.  In  terms  ot 
the  parameter  K,  a measure  of  the  rate  of  change  of  external  conditions  is 
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T^  e ds  Ot 


where  is  the  time  scale  over  which  eoge  conditions  change.  If  is  taken  as  re- 


presentative of  the  time  required  for  a boundary  layer  to  react  to  changing  condi- 
tions (i.e.,  to  achieve  equilibrium),  then  it  follows  that 


'e  " ^ 


for  d typical  nonequi  1 ibriuin  la.iiinarizing  flow.  This  observation  is  derived  from  the 
qualitative  features  of  the  data  and  is  also  the  basic  hypothesis  of  the  rec°nt  theory 
of  Deissler  (Reference  6). 


It  is  obvious  that  is  a complex  function  of  the  turbulent  velocity  field  in 


a boundary  laye*'.  In  fact,  it  represents  the  sum  total  of  the  dynamical  self- 
interactions between  turbulent  components  as  well  as  coupling  with  the  mean  velocity 


field.  A theoretical  prediction  of  would  require  the  solution  of  a complex  set  of 


partial  differential  equations,  which  is  a task  that  is  the  subject  of  current  basic 
research  and  obviously  beyond  the  scope  of  design  calculations 


Nevertheless,  correlation  of  the  low  speed  data  must  in  some  way  account  for 
nonequilibrium  effects.  This  is  shown  hj  an  application  of  classical  correlations 
to  the  experimental  conditions.  For  these  test  cases,  the  correlations  given  by 
White  (Reference  15)  were  chosen,  which  are 
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(log  Re  1 
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P = 0.8(0. 5 + 6) 


2 X 3.179P  + 1.5P^ 
''  ■ ' 0.41(1  + P) 


H = 


/2/C^  - A 


These  correlations  are  the  result  of  a large  amount  of  carefully  scrutinizeo 
daid  as  well  as  theoretical  considerations.  They  represent  a composite  of  the  work 
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of  Ludwig  and  Tillmanri  (Reference  16),  Coles  (Reference  17)  and  others.  IJasic  ele- 
ments in  the  construction  of  these  expressions  are  the  "law  of  the  wall"  and  the  "law 
of  the  wake",  whicti  are  described  by  White  (Reference  15)  and  > iJely  known.  It  should 
be  no"ed  that  these  are  usually  regarded  as  valid  for  equilibrium  boundary  layers  only, 


and  are  sometimes  referred  to  as  "equilibrium  correlations".  As  such,  the  usual  ap- 
plication of  these  equations  never  encounters  the  circumstances  where 


which,  for  the  large  K of  a laminarizing  flow,  is  very  common  and  results  in  ai  imagi- 
nary value  of  P.  Taking  into  account  the  fact  that  the  parameter  P has  its  origin  in 
the  "law  of  the  wake",  P has  been  set  equal  to  zero  for  any  1/2,  consistent  with 
the  theory  of  Oeissler  (Reference  6)  and  with  the  inference  from  experimental  velocity 
profiles  that  the  wake  component  disappears  with  large  pressure  gradient. 

lyi)i(ul  results  of  integration  of  the  integral  momentum  equation  using  these 
correlations  are  shown  in  Figure  7.  The  initial  conditions  for  the  integrations 
(Re,  (o))  were  taken  from  the  data,  as  were  the  distributions  of  edge  velocity.  Also 
shown  on  Figure  7 are  the  data,  which  are  seen  to  compare  very  poorly  with  the 
predictions.  The  data  show  a streaniwise  development  of  Re||  different  from  the  pre- 
dictions and,  more  importantly,  a very  different  trend  for  C^(s).  If  a boundary  layer 
is  considered  to  have  laminarized  at  that  strpaf,«ise  position  where  Re^,  and  attain 
values  equal  to  those  given  by  a laminar  prediction  (no  other  definition  is  possible 
with  integral  parameters),  then  Figure  7 shows  that  the  equilibrium  correlations  fail 
to  predict  "laminarization".  Further  examination  of  the  data  shows  that  the  stream- 
wise  position  where  laminarization  occurs  is  downstream  of  the  point  where  K is  a max- 
imum This  is  a manifestation  of  nonequilibrium  effects,  and  can  be  characterized  as 
a "delay". 

In  qualitative  terms,  some  of  the  data  show  a very  pronounced  delay  and  some 
a,'  not  while  other  data  indicate  a possible  approach  to  turbulent  sink  flow  behavior 
which  the  data  of  Blackwelder  and  Kovasznay  (Reference  8)  lack.  Therefore,  it  is 
necessary  to  identify  the  reasons  for  these  apparently  different  modes  of  laminari- 
zation and  to  identify  a parameter  that  correlates  these  effects. 

As  a preliminary  part,  of  this  task  and  as  a more  or  less  routine  check  on  the 
data  and  the  code  used  f^r  its  analysis,  the  measured  distributions  of  Re,  H,  K and 
Cj,  were  examined  for  consistency  with  respect  to  the  integral  momentum  equation.  The 
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Iiiethud  of  iiSllVJ  SiiiuuthcJ  CuLilC  Splllit  flt^  anj  Uunyo-KuUa  liilfcy.-allull , proviulisly 
describetl,  was  used.  The  results  are,  that  generally,  the  measured  and  calculated 
distributions  of  Re  ^ have  the  sasie  trends  but  differ  significantly  in  magnitude  as 
H isH'il  i_  U.o  id  sV.J  <:  ts.  5*1 

calculated,  respectively,  a convenient  measure  of  the  discrepancy  is  the  quantity 


(Re  - Re  )/Re 

' om  oc"  ijin 


An  annotated  summary  of  the  discrepancies  in  terms  of  the  above  parameter  for  each 
data  set  is  as  follows: 

A.  Blackwelder,  et  al.  (Reference  8) 

• 10  percent  discrepancy  in  region  of  niaxiimiin  laminarization 

• 50  percent  discrepancy  in  retransition  region 

8.  Patel,  et  al.  (Reference  9) 

« -93  percent  discrepancy  in  region  of  maxiimim  laminarization 

( Noticeably  different  trend  in  approach  to  miininum  Re  ^ 

C.  Badri  Narayanan,  et  al.  (Reference  10) 

• -53  percent,  -60  percent,  -90  percent  (worst  case)  increasing  with 
initial  Re,, 

« Better  agreement  i found  when  smoothing  is  applied  to  u^  and  K is 
recalculated 

Additional  understanding  of  the  inconsistencies  in  the  data  comes  from  rever- 
sing the  procedure.  That  is,  compute  dRe^yOx  from  a smoothed  cubic  spline  fit  to 
the  measured  Re,,  distributions,  use  measured  H and  Re  distributions  and  calculate 
from  the  inte<|ral  momentum  equation.  This  procedure  gives  Cj.  from  0.1  to  12  times 
the  iiieasured  values.  In  some  case,  for  certain  regions,  calculated  is  negative, 
whiih  is  a physical  imposstbi  lity. 

A number  of  cross  checks  of  the  code  and  the  method  of  analysis  were  made.  It 
was  shown,  toi-  example,  tnat  smoothing  of  the  measured  edge  velocity  distribution-, 
results  in  rressure  gradient,  (K),  distributions  slightly  different  than  those  cal- 
culated 1)/  the  authors  of  the  experimental  studies.  However,  such  differences  are 
insufficient  tu  account  for  the  overall  lack  of  comparison.  Another  test  case  was 
the  calculation  of  a flat  plate  boundary  layer  and  a separate  test  on  the  integration 
subi'outine  that  compared  the  numerical  solution  to  a known  exact  solution.  In  all 
cases  the  code  and  the  numerical  methods  were  found  to  bo  correct. 
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Figure  9.  Calculated  and  measured  Ren  distributions  for  the  data  of 
Patel  and  Head  (Reference  5). 


Measured 


ated  and  measured  Rea  distributions  for  the  data  of 
Narayanan  and  Ramjee” (Reference  10). 


A further  inquiry  examined  the  question  of  whether  or  not  the  integral  momen- 
tum equation  is  valid  for  a strongly  accelerated  flow.  The  only  aspect  of  this  ques- 
tion that  needs  to  be  addressed  concerns  the  boundary  layer  equations  themselves, 
because  the  von  Kantian  equation  is  just  an  integrated  form  of  these  equations.  The 
usual  boundary  layer  approximation  of  the  Navier-Stokes  equations  neglects  a term  in 
the  streamwise  momentum  equation,  i.e.,  it  is  assumed  that 

r'u  , a*u 


where  n is  normal  to  the  wall. 

This  inequality  was  examined  for  the  case  of  the  present  circumstances  and 
found  to  he  satisfied;  the  neglected  term  is  indeed  several  orders  of  magnitude 
smaller.  No  evidence  supporting  the  need  for  the  full  Navier-Stokes  equations  has 
been  found. 

With  regard  to  the  experiments,  there  are  a number  of  well  known  grounds  on 
wnicM  to  question  the  validity  of  data  of  this  kind.  In  the  present  case,  it  is  ap- 
propriate to  note  these; 

a.  The  shape  of  the  mean  velocity  profile  is  more  than  usually  difficult  to 
measure  accurately 

b.  Wall  shear  stress  measurements  by  ordinarv  methods  are  questionable 

r.  The  mean  ilow  may  not  be  two-dimensional. 

It  is  not  clear  which  one  or  more  of  the  above  factors  is  most  relevant  to  each  of  the 
experiments  considered.  One  of  the  authors  (Blackwelder,  Reference  18)  has  reported 
that  there  were  initial  doubts  about  the  two-dimensionality  of  the  flo.v  lot  that  the 
apiMratus  was  corrected.  The  data  from  this  and  another  experiment  wore  tested 
against  the  theory  of  Deissler  (Reference  19)  and  it  was  reported  by  Deissler  (Ref- 
erence 19)  that  the  wall  shear  stress  measurements  of  Blackwelder  and  Kovasznay  (Ref- 
erence 0 did  not  compare  as  well  as  did  those  of  Patel  and  Head  (Reference  9). 

Overall  1 ■■  seems  (air  to  judge  that  there  are  serious  fundamental  questions 
at'oat  tne  Tata  ttiat  are  presently  unanswerable.  Under  these  conditions,  it  is  only 
possible  to  get  ser  i-qiiantitative  infonnation  about  the  trend  of  the  data  and  to  make 
as  mncti  use  ot  this  as  possible. 

3.?  LAITIf.ARIZAI  ION  ON  NOStflP  CONFIGURATIONS 

( rom  ihe  results  ot  the  previous  sections,  it  is  evident  that  a definitive 
Jclineation  ot  the  mectianism  of  nosetip  laminarization  is  not  yet  available.  How- 
ever, It  is  shown  below,  that  a combination  of  an  empirical  expression  and  some 


re.isonable  arguments  leads  to  a useful  correlation  for  the  prediction  of  nosetip 
laminari^ation  in  design  studies.  The  arguments  are  partly  intuitive  but  also  based 
on  known  theoretical  and  experimental  results  for  flows  with  normal  acceleration 
(streamwise  curvature).  This  correlation  is  described  in  the  following,  starting 
with  its  empirical  basis  and  showing  the  results  of  test  cases  at  the  end  of  the 
section. 


3.2.1  The  Correlation 

The  degree  of  laminarization  at  the  tangent  point  of  a sphere-cone  was  shown 
by  Anderson  and  Jackson  {Reference  20)  to  be  well  c'rrelateu  by  conditions  at  the 
sonic  point.  When  the  degree  of  laminarization  is  defined  as 


where  g is  heat  flux  and  the  subscripts  T and  L refer  to  turbulent  and  laminar 
respectively,  it  v/as  ..hown  that 


fore  cone 


where 


in  which  T^/T^  is  the  ..'all  temperature  ratio,  u is  the  streamwise  velocity,  k is  the 
surface  rough. less  height,  is  the  streamwise  radius  of  curvature,  j is  the  kine- 
matic viscosity,  e refers  to  the  edge  of  the  boundary  layer,  and  * refers  to  the 
sonic  point.  The  correlation  and  the  data  are  shown  on  Figure  11. 

Ihe  data  base  for  tins  correlation  includes  data  from  the  rough  wall  calorim- 
eter experiments  of  the  PANT  program  (Series  A and  J)  and  those  from  the  NOL  ART 
tests  (Series  3 and  4).  These  data  cover  a fairly  wide  range  of  nose  radius,  sur- 
face roughness,  unit  Reynolds  number,  wall  temperature  ratio  and  geometry. 

The  needs  of  the  current  study  are  met  by  establishing  a procedure  for  cal- 
culating the  distribution  of  T from  an  appropriate  point  upstream  to  the  cone-tangent 
point  where  j.  is  known.  There  are  two  main  considerations  to  be  taken  into 
account  in  this  process.  The  first  is  the  choice  of  variables  and  the  second  is  the 
determination  of  an  analytical  form. 
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With  tegdi'c.!  to  the  choice  of  variables,  it  is  noted  that  on  a nosetip,  the 
pressure  gradient  parameter 


■'>-  dUg 


K = 


u„ 


is  infinite  at  the  forward  stagnation  point  and  decreases  iiionotonically  with  s, 
except  for  a region  somewhat  downstream  of  the  sonic  point,  where  it  reaches  a small 
relative  inaxiinum.  A typical  distribution  is  shown  in  Fig.  12,  along  with  the  cal- 
culated laminar  and  turbulent  and  the  measured  heat  transfer  distributions  for  one 
of  the  PANT  Series  J tests  (Run  602). 

The  point  to  be  noted  is  that  in  the  region  of  laminarization  on  a nosetip, 

K is  always  much  less  (by  a factor  of  approximately  10)  than  it  is  at  the  location 
of  transition.  In  contrast,  the  plane  flow  data  (with  streamwise  pressure  gradient) 
of  Hlackweldei'  and  Kovasznay  (Reference  8)  and  others  show,  at  least  qualitatively, 
that  laminarization  downstream  of  a point  where  a given  turbulent  condition  exists 
is  the  result  of  an  increase  of  K with  streamwise  distance.  Thus,  it  seems  unlikely 
that  nosetip  laminarization  is  strongly  influenced  by  the  streamwise  pressure  gra- 
dient. 

The  other  aspect  of  the  comparison  of  the  effect  of  streamwise  pressure  gra- 
dient in  low  speed  flows  with  the  case  of  a nosetip  is  the  magnitude  of  K.  Figure 
12  is  representative  of  the  circumstances  on  a nosetip,  and  shows  that  K • 3(10’’) 
over  most  of  the  sphere  (S/R|,j  < 1.2).  This  is  large  enough  to  cause  laminarization 
according  to  the  usual  low  speed  criterion  for  smooth  walls.  However,  the  effect  of 
roughness  on  low  speed  laminarization  is  not  known  (either  experimentally  or  theo- 
retically). It  is  intuitively  reasonable  that  the  value  of  K necessary  for  laini- 
narization  increases  with  increasing  roughness.  Hence,  it  is  argued  that  roughness 
is  a ma.ior  factor  in  determining  the  state  of  a boundary  layer  on  a nosetip  and 
accounts  for  the  differences  between  the  magnitude  of  K on  a nosetip  and  in  a low 
speed  flow  at  points  where  the  degree  of  laminarization  is  similar. 

The  remaining  variable  in  the  fore-cone  correlation  is  the  streamwise  radius 
nt  curvature,  Rj- , and  it  follows  that  the  extension  of  Ip  p should  involve  a varia- 
tilc  associated  with  this  effect.  The  work  of  Prandtl  (Reference  21),  Rotta  (Ref- 
erence 22),  .and  Thomann  (Reference  23)  supports  the  contention  that  the  main  effect 
of  streamwise  curvature  is  to  establish  a pressure  gradient  normal  to  the  surface, 
winch  inhibits  turbulent  mixing.  The  magnitude  of  the  normal  pressure  gradient  is 
conveniently  measured  by 


f>  u 

e e 


which  typically  has  a distribution  on  a nosetip  as  shown  in  Figure  12. 
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Since,  for  a plane  flow  with  pressure  gradient,  there  is  unequivocal  evidence 
for  a "lag"  effect,  and  s’  .,.e  the  data  of  Thoniann  (Reference  23)  suggest  a similar 
behavior  when  streamwise  curvature  is  present,  it  is  reasonable  to  conjecture  that 
some  aspect  of  this  "history"  of  the  normal  pressure  gradient  is  relevant.  A simple 
form  suggested  by  these  considerations  is 


Hs)  f 


where  S is  stroaimvise  distance  and  Sp  is  some  reference  location.  It  is  reasonable 
that  the  normal  pressure  gradient  measured  relative  to  that  at  the  reference  location 
should  be  in  the  integrand.  Thus,  it  is  postulated  that  the  integral 


. rVe’  l'W\  L 


is  proportional  to  the  degree  of  lannn:,rization  on  a nosetip.  Normalizing,  and 
imposing  the  requi r'en.ent  that  i'{Sp  c ) ' q Qives 

T.-  r US) 

V(S)  = 

WiUi  this  form  for  v,  it  is  still  necessary  to  choose  Sp  and  one  obvious 
choice  is 


" Transition  " ‘'u 


further  analysis  of  the  data  may  establish  a better  reference  location,  possibly 
the  sonic  point,  but  the  range  of  streamwise  position  in  which  Sp  may  be  taken  is 
not  extremely  lai-ge.  "^he  reason  for  this  restriction  is  that  it  appears  necessary 
that  I(S)  . 0 for  all  S ' Sp  , otherwise  it  is  doubtful  that  V(S)  will  follow  the 
trend  of  the  data. 


A sample  calculation,  with  the  code  modified  to  include  r(S),  is  shown  in 
figure  13  for  the  case  of  a PANT  Series  J test.  Figure  13  shows  reasonable  agree- 
;.ient  between  the  correlation  and  the  data,  as  do  other  test  cases.  The  only  obvious 
criticism  is  that  the  degree  of  laminarization  is  sligh’ly  overestimated  by  the  cor- 
relation in  the  region  where  S is  about  halfway  between  S^.  and  Sp  ^ . However,  it 
is  judged  that  I(S)  is  a reasonable  extrapolating  expression  and  useful  for  estimating 
nosetip  laminarization. 
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SECTION  4 
SENSITIVITY  STUDY 


In  order  to  assess  the  impact  of  relaminarization  on  reentry  nosetip  shape 
chantje,  the  model  described  in  the  previous  section  was  used  to  compute  histories  of 
heat  transfer  distributions  over  sphere  cone  nosetips  on  ballistic  trajectories.  Pre- 
dictions were  made  for  four  reentry  configurations,  with  and  without  relaminarization. 
Descriptions  of  th'  se  cases  are  given  in  Table  1.  All  cases  were  run  for  a peak  to 
valley  roughness  height  of  0.5  mils  which  is  nominally  characteristic  of  ATJ-S  graphite 
as  well  as  CE  2-2-3  carbon/carbon.  As  noted  in  Table  1,  the  same  typical  values  of 
tiic  ballistic  coefficient  and  entry  velocity  were  used  for  all  four  cases.  Since  nose 
radius  and  entry  angle  were  expected  to  influence  the  overall  impact  of  laminarization 
M'ost  directly,  predictions  were  made  for  two  values  of  each  of  these  parameters.  The 
ti  ll  jectorios  for  the  two  entry  angles  are  shown  in  Figure  14. 


Because  the  relaminarization  model  is  most  valid  for  sphere  cone  geometries 
and  because  this  geometry  is  quite  characteristic  of  the  actual  shape  for  the  portion 
of  tiie  trajectory  during  which  laminarization  is  most  important,  the  original  geo- 
metry was  assumed  for  the  entire  trajectory.  This  approximation  may  slightly  under- 
estimate the  influence  of  laminarization  because  the  corner  radius  decreases  in  the 
presence  of  turbulent-relaminarized  flow.  The  Sandia  Nosetip  Analysis  Procedure  (SNAP) 
code  was  run  without  ablation,  with  the  wall  temperature  assigned  at  7000"R. 


Transition  onset  altitude  and  subsequent  transition  locations  are  computed 
in  tiie  SNAP  code  using  the  rough  wall  boundary  layer  transition  criteria  developed 
under  Uie  PANT  program  (see  Reference  24). 


255*,  onset 
215,  location 
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and 
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Evaluated  at  the  sonit.  i;oint 
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TABLE  1 


BALLISTIC  AND  GEOMETRIC  PARAMETERS 
USED  IN  SENSITIVITY  STUDY 


case 

1' 

Ibf/ff’ 

Ve 

ft/sec 

>e 

degree 

Rn 

inch 

Trans i Lion 
A1 ti tude 
ft 

1 

3700 

22,000 

-26 

0.75 

58,000 

2 

3/00 

22,000 

-26 

1.50 

59,000 

3 

3700 

22,000 

-40 

0.75 

60,000 

4 

3700 

22,000 

-40 

1.50 

61 ,000 
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Since  ablation  was  not  included  in  the  present  calculations,  ii'  = 0 and  i for  pure 
air  were  used  in  coiiiputing  iji.  This  approximation  yields  values  of  i/  which  are  10  per- 
cent low  compared  to  values  computed  including  the  effects  of  ablation  (see  Figure  A-20, 
Reference  24).  This  in  turn  implies  that  the  transition  altitudes  reported  herein  are 
about  2000  f-  above  those  which  would  be  calculated  with  ablation.  This  has  an  in- 
significant effect  on  the  calculation  of  relaminarization  at  any  altitude  and  only  a 
small  effect  on  the  integrated  influence  of  relaminarization  over  the  reentry  tra- 
jectory. 

Relaminarization  has  its  greatest  influence  on  the  heat  transfer  distribution 
around  a nose',  i;..  at  altitudes  just  below  transition  onset.  Figure  15  shows  predicted 
heat  transfer  distributions  for  case  1 at  53,000  ft,  about  5,000  ft  below  transition 
altitude.  At  the  cone  tangent  point,  ip  = 0.77,  which  is  its  maximum  on  the  nosetip. 

Figures  16,  17,  18  and  19  show  predicted  heat  transfer  distributions  for  the 
four  cases,  all  at  the  same  altitude,  41,000  ft.  As  expected  frimi  the  cone  tangent 
point  I'elaminarization  correlation,  'ncreasing  the  nose  radius  decreases  the  in- 
fluence of  relaminarization.  It  is  also  apparent  from  these  four  fit|ures  that  for 
the  same  ballistic  coefficient,  entry  velocity  and  nose  radius,  the  entry  angle  has 
only  a minor  effect  upon  laminarization  at  this  altitude  where  the  velocities  on  the 
two  trajectories  are  nearly  the  same. 

Figures  20,  21,  22  and  23  indicate  how  the  predicted  influence  of  laminariza- 
tion on  the  cone  tangent  point  heat  transfer  varies  with  altitude  for  the  four  cases 
studied.  The  nose  radius  eff»ct  is  as  expected.  Although  the  heat  transfer  on  the 
steeper  trajectory  ■=  -40  ) is  considerably  higher,  the  fractional  reduction  in 
neat  transfer  associated  with  relaminarization  at  any  altitude  is  about  the  same  for 
the  two  trajectories. 

Figures  24  through  27  show  the  predicted  influence  of  laminar! zation  on  the 
total  fieat  transfer  at  the  cone  tangent  point.  Focusing  on  an  altitude  of  30,000 
ft,  cases  one  tfirough  four  show  ' reduction  of  integrated  heat  transfer  at  the  cone 
tangent  point  of  27,  17,  28  and  17  percent,  respectively. 

The  reductions  in  heat  transfer  will  be  directly  reflected  in  reductions  in 
recession.  Although  fore  cone  ablation  is  not  crucial  from  a recession  point  of  view, 
this  region  of  the  nosetip  is  critical  to  thermostructura  1 performanci' . The  effects 
of  laminarization  on  heat  transfer  are  probably  important  with  respect  to  thermal 
stress  considerations  and  an  assessment  of  its  thermnstructural  impact,  which  was  not 
within  the  scope  of  the  (iresent  work,  should  be  carried  out. 
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Figure  21.  Fore  cone  heat  transfer  rate  as  a function  of  altitude 
Case  2. 


Turbuleri 


SECTION  5 


CONCLUSIONS 


Tiil'  >;nrn."ltiti()it  '.tudy  of  relaininarization  of  piane  flows  due  lo  slrear.wtse 
pressiii'e  i|r.ulU'Ml  his  ‘howii  the  existence  of  inconsistencies  in  the  experiinenldl  data 
Hence,  a full  account  of  the  effects  of  streaniwiso  acceleration  and  non-e(]ui  1 ibrium  are 
not  possible  without  further  theoretical  and  experimental  work  to  resolve  these  in- 
consistencies. However,  qualitative  and  ' -quantitative  arguments  have  been  pre- 
sented that  indicate  that  the  iiiechanisin  of  ..usetip  laininarization  is  most  likely 
associated  with  the  effects  of  surface  roughness  and  normal  acceleration.  As  a 
corollary  to  this  result,  it  has  been  shown  that  the  commonly  used  streamwise  p'  js- 
sure  gradient  criterion  for  laminarization  does  not  give  a good  correlation  of  wind 
tunnel  data  on  nosetip  configurations. 


An  extension  of  an  existing  correlation  for  the  degree  of  laminarization  at 
tlie  sphere-cone  tangent  point  of  a nosetip  has  been  developed.  This  extended  cor- 
relation predicts  the  distribtuion  of  heat  transfer  rate  from  the  point  of  transi- 
tion to  downstream  of  the  sphere-cone  tangent  point  by  means  of  an  integral  of  a 
measure  of  the  normal  pressure  gradient. 


Heat  transfer  distributions  on  reent’"y  nosetips  with  two  nose  radii  and  on 
two  trajectories  have  been  calculated  in  a sensitivity  study  of  the  effect  of  re- 
laminarization.  It  has  been  shown  that  the  time  integral  of  the  heat  transfer  rate 
on  the  fore  cone  o*’  a nosetip  is  significantly  decreased  by  laminarization.  The  re- 
duced local  heat  transfer  rates  that  result  from  laminarization  are  of  importance 
primarily  to  the  thermoslructural  performance  of  a nosetip,  and  the  present  results 
indicate  the  need  to  make  a detailed  study  of  the  thermostructural  ramifications  of 


these  effects. 
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